Energy and expectation values of the PsH system. 
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Close to converged energies and expectation values for PsH are computed using a ground state 
wave function consisting of 1800 explicitly correlated gaussians. The best estimate of the PsH°° 
energy was -0.789196740 hartree which is the lowest variational energy to date. The 27 annihilation 
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PACS numbers: 36.10.Dr, 36.10.-k, 31.15.Ar 

The calculation of positronium hydride (PsH) repre- 
sents one of the simplest possibilities for studying mixed 
electronic and positronic systems. Since its stability was 
first identified in 1951 by Ore a variety of methods 
have been applied to determine its structure. These in- 
clude variational calculations with Hylleraas type basis 
sets 0, 0, 0, H|, variational calculations with explicitly 
correlated gaussians (ECGs) El E I 11101, quantum Monte 
Carlo methods jlTjL fill IT2I Il3lll4| and most recently the 
configuration interaction method 

lowest variational energy for PsH°° prior to the present 
article was that of Yan and Ho 3 • Their largest calcula- 
tion gave an energy of -0.7891967051 hartree. Bubin and 
Adamowicz used a 3200 dimension ECG basis to give an 
energy of -0.788870707 hartree for PsH 1 |l|. 

In this work, the stochastic variational method (SVM) 
is used to construct a wave function with a lower energy 
that the best wave function of Yan and Ho. Indeed, the 
best SVM energy of -0.789196740 hartree is even lower 
than the value estimated by Yan and Ho as the varia- 
tional limit (e.g. -0.7891967147(42) hartree). 

The SVM used for this work has been described in 
a number of articles [20I I2H 122^ and only the briefest 
description is given here. The SVM expands the wave 
function in a linear combination of ECGs. Such basis 
functions have Hamiltonian matrix elements that can be 
computed very quickly and the energy is optimized by 
performing a trial and error search over the exponential 
parameters that define the basis. The SVM has been used 
to solve a number of many-body problems in different 
areas of physics [3 H3 • 

For the present set of calculations a basis containing 
1800 ECGs was used for the final calculation. All the 
optimizations of the ECG basis were done with the H 
mass set to 00. The annihilation rates given in Tables [I] 
and [H] are proportional to the probability of finding an 
electron and a positron at the same position in a spin 
singlet state according to 

r = 47rr&*l£0&*(r< (1) 

i 

= 1.009394 x 10 11 ^2(S(r t - r p )) s , (2) 

i 

0,13, 13] • The sum is over the electron coordinates, the 
8- function expectation is evaluated in afj, and V is given 



numerically in s _1 . The operator Of p is a spin projection 
operator to select spin singlet states for the ip electron- 
positron pair. 

Table [I] lists a number of expectation values obtained 
from a sequence of increasingly larger calculations. The 
net energy improvement when the basis was increased 
from 900 to 1800 ECGS, while being subjected to addi- 
tional optimization, was 1.98xl0~ 7 hartree. It is worth 
noting that the energy of the largest calculation, namely 
-0.789196740 hartree, is lower than the previous best en- 
ergy of Yan and Ho 0, namely -0.7891967051 hartree. 
Yan and Ho examined the convergence pattern associ- 
ated with their sequence of increasingly larger calcula- 
tions and estimated that the true energy was actually 
9.6(4.2) x lO" 9 hartree lower (e.g. -0.7891967147(42) 
hartree). The present calculation indicates that the ac- 
tual correction should have been more than three times 
as large as that estimated by Yan and Ho. Although the 
sign of size of energy correction is not large, it is apparent 
that the procedure used to determine the energy correc- 
tion is faulty. In Hylleraas calculations one typically does 
some sort of non-linear optimization to choose the expo- 
nential parameters that give the minimum energy. This 
has the unintended byproduct of distorting the conver- 
gence pattern of the energy and thus introdu cing large 
uncertainties in the extrapolation of the energy [25] . This 
problem is probably more widespread than just the PsH 
calculation of Yan and Ho. It could occur whenever one 
extrapolates a sequence of energies while using a family 
of basis functions that are characterized by a parameter 
which has been subjected to a non-linear optimization. 

The coalescence matrix elements, (S(e~ — e~)) and 
(<5(H + — e + )} were more sensitive to the increase in ba- 
sis size than any other quantity. This sensitivity is due 
to the fact that the wave function amplitude between 
two repelling particles is expected to be small at their 
coalescence point and the ECG functional form is not 
the natural choice to describe the behavior of the rel- 
ative wave function for two strongly repelling particles. 
With respect to the more physically interesting observ- 
ables, the annihilation rate, T varied most as the basis 
dimension was increased. But, the increase in T was just 
larger than 0.1% when the basis was increased from 900 
to 1800. 

A comprehensive set of the best set of expectation val- 
ues are listed in Table lH| They are compared with the 
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TABLE I: Behavior of some PsH expectation values for a sequence of ECG type variational calculations of increasing size. 
All quantities are given in atomic units with the exception of the 27 annihilation rates which are in units of 10 9 s _1 . Some of 
the data for the earlier calculation |7J have not been published before, the data attributed to these calculations were computed 
using the same ECG basis. 
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(V)I(T) + 2 



E 



750 [7J 3.661596 7.812895 0.3705556 

900 3.661613 7.812961 0.3705554 

1200 3.661621 7.813024 0.3705550 

1500 3.661624 7.813040 0.3705549 

1800 3.661624 7.813046 0.3705549 
Hylleraas N = 5741 [4, 5] 
Hylleraas N — » 00 extrapolation [4, 5] 



3.480249 4.39845 x 10~ 3 1.63863 x 10" 

3.480263 4.39321 x 10" 3 1.63635 x 10" 

3.480270 4.38188 x 10~ 3 1.63153 x 10" 

3.480271 4.37628 x 10~ 3 1.62850 x 10" 

3.480272 4.37639 x 10 -3 1.62828 x 10" 



2.46852 
2.46879 
2.47129 
2.47134 
2.47178 
2.47258 
2.47264(2) 



5.51 x 10"' 
7.96 x 10~ 7 
2.21 x 10~ 7 
1.30 x 10~ 7 
7.3 x 10" 8 



-0.789195993 
-0.789196542 
-0.789196673 
-0.789196718 
-0.789196740 
-0.789196705 
-0.789196715(5) 



results of another, but completely independent, large ba- 
sis SVM calculation 0. The expectation value for the 
virial theorem (V) / (T) provides an estimate of the wave 
function accuracy and the deviation of (V) / (T) from -2 
was only 7.3 x 10~ 8 hartree. 

The energies of the different mass variants of PsH were 
computed by rediagonalizing the Hamiltonian with the 
same basis but with m H i set to 1836.1527 m e , m^-i set 
to 3670.483 m e and m H 3 set to 5496.899 m e . The ener- 
gies of PsH 1 , PsH 2 and PsH 3 were -0.788870618 hartree, 
-0.789033556 hartree and -0.789087767 hartree respec- 
tively. The energy of the 3200 ECG wave function of 
Bubin and Adamowicz 19] for PsH 1 was -0.788870707 
hartree, which is 1.0 x 10 -7 hartree below the present 
energy. 

To summarize, a close to converged binding energy is 



reported for the PsH°° ground state. The present energy 
is 2.5 x 10 -8 hartree lower than the estimated variational 
limit of Yan and Ho. The procedure by Yan and Ho to 
estimate the variational limit probably tends to underes- 
timate the size of the necessary energy correction. 

Although the present energy is better than that of Yan 
and Ho, this does not necessarily mean that the present 
SVM annihilation rate is more accurate. Any basis of 
ECGs (which cannot satisfy the exact inter-particle cusp 
conditions) will have a tendency to underestimate the 
electron-positron coalescence matrix element. Table [I] 
shows a consistent increase in T as the size of the cal- 
culation in increased. 
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TABLE II: Properties of the PsH ground state. Data are 
given for H assuming infinite mass. All quantities are given 
in atomic units with the exception of the annihilation rates 
which are in units of 10 9 s _1 . The positron and electron 
kinetic energy operators are written as T+ and T_. 



Property 
N 

(V)/(T) + 2 
E 

(T-) 
(T+) 

( r H+e-) 
(»H+e+) 

v%- e -) 
(r e + e -) 

(l/r-H+e-) 

(l/r H +e+) 

(lAe+e-> 

(r 2 e - e -) 
(r 2 +e ~) 

(1Ah +e+ ) 

(lAe- e+ > 

(*(H+-e-)) 
(<5(H+-e+)> 
(5(e"- e -)) 
(5(e+-e-)) 



Present SVM 
1800 
7.3 x 10" 8 
-0.789196740 
0.3261733 
0.1368503 
2.311526 
3.661624 
3.574787 
3.480272 
0.7297090 
0.3474618 
0.3705549 
0.4184961 
7.813046 
16.25453 
15.87546 
15.58427 
1.207067 
0.1721631 
0.2139099 
0.3491440 
0.177279 
1.62828 x 10~ 3 
4.37639xl0~ 2 
0.0244877 
2.47178 



SVM [9] 

1600 
6 x 10 -7 

-0.789165554 
0.3261732 
0.1368501 
2.311525 
3.661622 
3.574783 
3.480271 
0.7297087 
0.3474618 
0.3705549 
0.4184960 
7.813015 
16.25448 
15.87538 
15.58423 
1.207063 
0.1721637 
0.2139106 
0.3491428 
0.177186 
1.63857 x 10" 3 

4.3867xl0 -3 
0.024461 
2.46909 



